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Abstract 

We study stochastic tree fluid networks driven by a multidimensional Levy process. We 
are interested in (the joint distribution of) the steady-state content in each of the buffers, 
the busy periods, and the idle periods. To investigate these fluid networks, we relate the 
above three quantities to fluctuations of the input Levy process by solving a multidimensional 
Skorokhod reflection problem. This leads to the analysis of the distribution of the componentwise 
maximums, the corresponding epochs at which they are attained, and the beginning of the first 
last-passage excursion. Using the notion of splitting times, we are able to find their Laplace 
transforms. It turns out that, if the components of the Levy process are ‘ordered’, the Laplace 
transform has a so-called quasi-product form. 

The theory is illustrated by working out special cases, such as tandem networks and priority 
queues. 
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1 Introduction. 

Prompted by a series of papers by Kella and Whitt [20l l22l 25 , 23, there has been a considerable 
interest in multidimensional generalizations of the classical storage model with nondecreasing Levy 
input and constant release rate EH Ch. 4]. In the resulting networks, often called stochastic fluid 
networks, the input into the buffers is governed by a multidimensional Levy process. Recently, 
motivated by work of Harrison and Williams on diffusion approximations mm, the presence 
of product forms has been investigated mi maun mm. Recall that the stationary buffer-content 
vector has a product form if it has independent components, meaning that the distribution of this 
vector is a product of the marginal distributions. 

The results in these papers show that, apart from trivial cases, the stationary buffer-content 
vector of stochastic fluid networks never has a product form. Despite this ‘negative’ result, we 
show that it may still be possible to express the joint distribution of the buffer content in terms 
of the marginal distributions. This fact is best visible in the Laplace domain. For certain tandem 
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queues, for instance, the Laplace transform is a product that cannot be ‘separated’; we then say 
that the buffer-content vector has a quasi-product form. 

In the literature on stochastic fluid networks, there has been a focus on the stationary buffer- 
content vector W or one of its components. Here, we are also interested in the stationary distri¬ 
bution of vector of ages of the busy periods B and idle periods I. The age of a busy (or idle) 
period is the amount of time that the buffer content has been positive (or zero) without being zero 
(positive). Knowing these, it is also possible to find the distribution of the remaining length of the 
busy (or idle) period and the total length of these periods. 

We are interested in W, B, and I for a class of Levy-driven fluid networks with a tree structure, 
which we therefore call tree fluid networks. Our analysis of these networks relies on a detailed 
study of a related multidimensional Skorokhod reflection problem (see, e.g., Robert ,3S3)- Using its 
explicit solution, we relate the triplet of vectors ( W. B,I ) to the fluctuations of a multidimensional 
Levy process X. We also prove that the stationary distribution of the buffer-content vector is 
unique. 

Since our analysis of fluid tree networks is based on fluctuations of the process X, this paper 
also contributes to fluctuation theory for multidimensional Levy processes. Supposing that each 
of the components of X drifts to —oo, we write X for the (vector of) componentwise maximums 
of X, G for the corresponding epochs at which they are attained, and H for the beginning of the 
first last-passage excursion. Under a certain independence assumption, if the components of G are 
‘ordered’, we express the Laplace transform of ( X,G ) in terms of the transforms of the marginals 
( Xj , Gj). Since Xj is a real-valued Levy process, the Laplace transform of ( Xj , Gj) is known if Xj 
has one-sided jumps; see for instance Bertoin JU Thm. VII.4], 

We also examine the distribution of H under the measure Pj|,, which is the law of X given that 
the process Xf ; stays nonpositive. There exists a vast body of literature on (one-dimensional) Levy 
processes conditioned to stay nonpositive (or nonnegative), see the recent paper by Chaumont and 
Doney _SJ for references. Under the measure Pi, we also find the transform of ( X,G ). As a special 
case, we establish the Laplace transform of the maximum of a Levy process conditioned to stay 
below a subordinator, such as a (deterministic) positive-drift process. 

By exploiting the solution of the aforementioned Skorokhod problem, the results that we obtain 
for the process X can be cast immediately into the fluid-network setting. For instance, the knowl¬ 
edge of (X, G ) allows us to derive the Laplace transform of the stationary distribution of (W, B) in 
a tandem network and a priority system if there are only positive jumps, allowing Brownian input 
at the ‘root’ station. That is, we characterize the joint law of the buffer-content vector and the 
busy-period vector. With the Pi-distribution of H, we establish the transform of the idle-period 
vector I for a special tandem network. Our formulas generalize all explicit results for tandem fluid 
networks that are known to date (in the Laplace domain), such as those obtained by Kella [201 an d 
more recently by Dqbicki, Mandjes and van Uitert [ 0 ]. Most notably, quasi-products appear in our 
formulas, even for idle periods. 

To derive our results, we make use of the notion of splitting times. These essentially allow us 
to reduce the problem to the one-dimensional case. For real-valued Markov processes, splitting 
times have been introduced by Jacobsen CHI- Splitting times decompose (‘split’) a sample path of 
a Markov process into two independent pieces. A full description of the process before and after the 
splitting time can be given. However, since the splitting time is not necessarily a stopping time, 
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the law of the second piece may differ from the original law of the Markov process. We refer to 
Millar |,'iH IHS ; for further details, and to Kersting and Memi§oglu m for a recent contribution. 

The idea to use splitting times in the context of stochastic networks is novel. The known 
results to date are obtained with Ito’s formula m , a closely related martingale m, or differential 
equations ESI- Intuitively, these approaches all exploit a certain harmonicity. However, the results 
of Kyprianou and Palmowski m already indicate that there is a relation between these approaches 
and splitting. Splitting has the advantage that it is insightful and that proofs are short. Moreover, 
it can also be used for studying more complicated systems |E|. 

This paper is essentially divided into two parts. In the first part, consisting of Sections El El 
we analyze the fluctuations of an n-dimensional Levy processes X. The notion of splitting times 
is formalized in Section [3 These splitting times are first used to study the distribution of (X,G) 
in Sectional and then to analyze the distribution of H under Pi in Section 0] The second part of 
this paper deals with fluid networks. Section |B1 ties these networks to fluctuations of X , so that the 
theory of the first part can be applied in Sectional Finally, in Appendix^ we derive some results 
for compound Poisson processes with negative drift. They are used in Section El 

2 Splitting times. 

This paper relies on the application of splitting times to a multidimensional Levy process. After 
splitting times have been introduced, we study splitting at the maximum ('Section 12.1 1) and splitting 
at a last-passage excursion (Section 12.21) . 

Throughout, let X = (Ad,..., X n )' be an n-dimensional Levy process, that is, a cadlag process 
with stationary, independent increments such that A(0) = 0 € R n . Without loss of generality, as in 
Bertoin [Ij, we work with the canonical measurable space (12, JF) = (D([0, oo), M d U {9}), £>), where 
B is the Borel er-field generated by the Skorokhod topology, and d is an isolated point that serves 
as a cemetery state. In particular, X is the coordinate process. Unless otherwise stated, ‘almost 
surely’ refers to P. All vectors are column vectors. 

The following assumption is used extensively throughout this paper: 

D X k (t) —> —oo almost surely, for every k. 

We emphasize that a dependence between components is allowed. In the sequel, X k (t) (or 
X k (t)) is shorthand for sup s<t X k (s) (or inf s < t X k (s)). Due to D, X k := X k (oo) is well-defined 
and almost surely finite for every k. Furthermore, we write X = (X \,..., X n )'. 

The following two definitions are key to further analysis. The second definition is closely related 
to the first, but somewhat more care is needed on a technical level. Intuitively, for the purposes of 
this paper, there is no need to distinguish the two definitions. 

Definition 2.1 We say that a random time T is a splitting time for X under P if the two processes 
{A(t) : 0 < t < T} and {X(T + t) — X(T) : t > 0} are independent under P. We say that 
T is a splitting time from the left for X under P if the two processes {X(t) : 0 < t < T} and 
{X{T + t) — X (T—) : t > 0} are independent under P. 

Note that if X is a Levy process under P with respect to some filtration X which includes 
the natural filtration, any JF-stopping time t is a splitting time for X under P. In fact, the Levy 
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assumption implies that {X(t + t) — X(t) : t > 0} is not only independent of {X(t) : 0 < t < t}, 
but that it also has the same distribution as {X(t) : t > 0}. 

We need some notions related to the initial behavior of X. For k = 1 ,n, set R k = inf{t > 
0 : X k {t) = X k (t)}. Since {X k (t) — X k (t) : t > 0} is a Markov process under P with respect to 
the filtration generated by X (see Proposition VI. 1 of [5]), the Blumenthal zero-one law shows that 
either R k > 0 almost surely (0 is then called irregular {or {X k (t.) — X k (t) : t > 0}) or R k = 0 almost 
surely (0 is then called regular for {X k (t) — X k {t) : t > 0}). We also set R k = inf{f > 0 : X k (t) = 
X fc (f)}, and define regularity of 0 for fX k (t) — X k (t) : t > 0} similarly as for {X k (t)—X k (t) : t > 0}. 
If R k = 0 almost surely, we introduce 

S k = S k := inf {t > 0 : X k {t) + X k (t)}. 

Again, either S k = 0 almost surely (0 is then called an instantaneous point for {X k (t) — X k (t) : 
t > 0}) or S k > 0 almost surely (0 is then called a holding point for {X k (t) — X k (t) : t > 0}). One 
defines S_ k , instantaneous points, and holding points for {X k (t) — X_ k (t) : t > 0} similarly if R k = 0. 


2.1 Splitting at the maximum under P. 

Let G k = G k := inf(t > 0 : X k (t) = X k or X k (t—) = X /,.} be the (first) epoch that X k ‘attains’ its 
maximum, and write G = (Gi,..., G n ) . Observe that G k is well-defined and almost surely finite 
for every fc by D. 


Lemma 2.1 Consider a Levy process X that satisfies D. 

(i) If Rk > 0 P -almost surely or X k is a compound Poisson process, then G k is a splitting time 
for X under P. 

(ii) If R k = 0 P -almost surely but X k is not a compound Poisson process, then G k is a splitting 
time from the left for X under P. 

Proof. We use ideas of Lemma VI.6 of Bertoin {3j, who proves the one-dimensional case under 
exponential killing. 

We start with the first case, in which the ascending ladder set is discrete. Set to = 0 and define 
the stopping times r n+ i = inf{t > r n : X k (t) > X k (t —)} for n > 0. Write N = sup{n : r n < oo}. 
Note that D implies that N < oo almost surely. 

Let F and K be bounded functionals, and apply the Markov property to see that for n £ Z + , 


E[F(V(f),0 < t < G k )K(X(G k +t) - X(G k ),t > 0 );N = n\ 


= E 


F(X(t),0 < t < T n )l {N > n} K(X(T n + t) - X(r n ),t > 0)l {sup Xk (t)=x k (Tn)} 


= E [F(A'(t),0 <t < r n )l{ 

AT>n}] ® 

= E [F(X(t),0 < t < r n )l {Ar > n }] E 


A (X(r n + t) X(r n ),t > 0)lj suPt ^^ x k (t)=x k (r„)} 
K (X(t),t> 0)l {supt x h (t)=o } ■ 


Summing over n shows that the processes { X(t ) : 0 < t < G k } and {X(G k + 1) — X(G k ) : t > 0} 
are independent. 

The argument in the case R k = 0 is more technical, but essentially the same. The idea is to 
discretize the ladder height structure, for which we use the local time J k at zero of the process 
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{X k (t) — X k (t) : t > 0}; see Bertoin jlj Ch. IV] for definitions. Note that i k ( oo) < oo almost surely 
by Assumption D. 

Therefore, we fix some e > 0, and denote the integer part of e _1 ffc(oo) by n = |_e _1 Ifc(oo)J. A 
variation of the argument for R k > 0 (using the additivity of the local time) shows that {X(t) : 
0 < t <I k ( ne )} and {X(£ k (ne) + t) — X(£ k (ne)) : t > 0} are independent. According to I3J 
Prop. IV.7(iii)], I k (ne) j G k as e j 0, which proves the lemma. □ 


2.2 Splitting at a last-passage excursion under P^. 

Let H k = H k := inf{f > 0 : sup s>t X k (s) / X k (t)} be the beginning of the first last-passage 
excursion, and write H = (Hi ,..., H n )'. 

In this subsection, we study the splitting properties of H k for some fixed k = l,...,n. We 
suppose that 0 is a holding point for {X k (t) — X_ k (t) : t > 0}, i.e., that R k = 0 and S_ k > 0 P-almost 
surely. Under this condition, the event {X k = 0} has strictly positive probability. Therefore, one 
can straightforwardly define the conditional law P[. of X given X k = 0. 

It is our aim to investigate splitting of H k under Pjj., but we only have knowledge of X under 
P. As a first step, it is therefore useful to give a sample path construction of the law Pf on the 
canonical measurable space (f2,.F). For this, we define a process X k 1 by 


X kl (t) 


X(t) 

X(R% ) )-X((R® + S$ ) 


t)~) 


if t G 
if t E 


p(i) oO)') . 
Rk ’Rk I i 

cU) pO)) 
Rk -—k ) ) 


(1) 


where R ^ 0) = 0, and for j > 1, 

:= inf {t > rig-V : X k (t) + X k (t)} , R^ := inf [t > : X k (t) = X k (t)} . 

In other words, X k 1 is constructed from the coordinate process X by ‘reverting’ the excursions of 
{X k (t)~X k (t):t> 0}. 

We have the following interesting lemma, which is the key to all results related to Pj;. For the 
random-walk analogue, refer to Doney jHj. 


Lemma 2.2 Consider a Levy process X that satisfies D. If R k = 0 and S_ k > 0 P -almost surely, 
then X k 1 has law P|, under P. 

Proof. Observe that R k > 0, and that the post-maximum process {X(G k + 1) — X(G k ) : t > 0} 
has distribution Pjj. (a proof of this uses similar arguments as in the proof of Lemma 12.11 see 
Millar (H E2j for more details). 

Fix some q > 0, and let e q be an exponentially distributed random variable, independent of X 
(obviously, one must then enlarge the probability space). The first step is to construct the law of 
{X(G g k + t) - X(G q k ) : 0 < t < e q - G q k }, where G q k := inf{f < e q : X k (t) = X k (e q ) or X k (t~) = 
Xk( e q)}- By the time-reversibility of X Lem. II.2], it is equivalent to construct the law of 
{X(F q ) - X((F q - t)-) : 0 < t < F q }, where F q := sup{t < e, : X k (t) = X k (e q ) or X k (t~) = 
X k (e q )}. 

To do so, we use ideas from Greenwood and Pitman E3j. Let £ k be the local time of {X k (t) — 
X_ k (t) : t > 0} at zero (since R k = 0, S k > 0, we refer to Bertoin [3J Sec. IV.5] for its construction). 
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Its right-continuous inverse is denoted by I k 1 . The X-excursion at local time s, denoted by X s , is 
the cadlag process defined by 

X s (u) := X {{e k \s-) + u ) A £\s)) - X s -)-) , u > 0. 

If Cf l (s —) = then we let X s be d , the zero function that serves as a cemetery. Since 

{X s : s > 0} is a cadlag-valued Poisson point process as a result of D, one can derive (e.g., with 
the arguments of Lemma II.2 and Lemma VI.2 of Hj) that the process 

W := (W(s) = {D(s),X s ) : s > 0} 

is time-reversible, where D(s ) := X (^ 1 (s)). After setting cr q := it can be seen that this 

implies that {(D(s), X s ) : 0 < s < o q } and {(D(cr q —) — D((cr q — s)—), X a i- S ) : 0 < s < a q } have the 
same distribution. In other words, one can construct the law of {X(F?)~ X((F k — t)— ) : 0 < t < F k } 
from the law of {X(t) : 0 < t < F k } by ‘reverting’ excursions as in 0- 

It remains to show that this construction is ‘consistent’ in the sense of Kolmogorov, so that one 
can let q —> 0 to obtain the claim. For this, note that the family {cr q } can be coupled with a single 
random variable through a q = l k 1 (ei/q). □ 

We now study the splitting properties of using the alternative construction of P|. given in 
Lemma 12.21 Since S_ k ^ is a P-stopping time with respect to the (completed) natural filtration of X, 
the Markov property of X under P with respect to this filtration [IJ Prop. 1.6] immediately yields 
the following analogue of Lemma , 12.11 

Lemma 2.3 Consider a Levy process X that satisfies D. If R k = 0 and S_ k > 0 P -almost surely, 
then H k is a splitting time for X under Pj^. Moreover, it has an exponential distribution under Pjj.. 

We remark that the construction and analysis of p! is the easiest under the assumption that 
R k = 0 and S_ k > 0 P-almost surely, which is exactly what we need in the remainder. A vast body 
of literature is devoted to the case n = 1, and the measure Pj is then studied under the assumption 
that R\ = 0. This is challenging from a theoretical point of view, since the condition that the 
process stays negative has P-probability zero. Therefore, much more technicalities are needed to 
treat this case. We refer to Bertoin .3 and Doney [2, for more details. See also Chaumont and 
Doney [5]. 

3 The P-distribution of (X,G). 

The aim of this section is to find the Laplace transform of the distribution of (X,G), assuming 
some additional structure on the process X. Thus, in the sequel we write X k -< Xj if there exists 
some K k j > 0 such that Xj — K k jX k is nondecreasing almost surely. 

Lemma 3.1 Suppose the Levy process X satisfies D. If X k -< Xj, then G k < Gj. 

Proof. First note that G k ,Gj < oo as a consequence of D. To prove the claim, let us assume 
instead that Gj < G k while X (t) := Xj(t ) — CX k (t) is nondecreasing for some arbitrary C > 0. 
Suppose that X k {Gk) = X k and Xj(Gj) = Xj ; the argument can be repeated if, for instance, 
Xk(G k —) = Xfc. The assumption Gj < G k implies that 

0 < X(G k ) - X(Gj) = Xj{G k ) -Xj-C [X k - X k (Gj)\ < 0, 
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meaning that X^ = X^(Gj). This contradicts Gj < G^ in view of the definition of G^. □ 

The following proposition expresses the distribution of ( X , G) in terms of those of ( X(Gk ), Gfc) 
and (X(Gfc—), Gfc). We denote the scalar product of x and y in M” by (x,y), and we write ‘cpd 
Ps’ for ‘compound Poisson’. Throughout this paper, the expression fl, a j x n/ x T should be 

read as (n, x (n. ; Pj) x 7 - 


Proposition 3.1 Suppose that X is an n-dimensional Levy process satisfying D and that X\ -< 
X -2 X n . Then for any a, (3 € M", 


Ee -(a,G)-(p,X) 


n —1 


n 


a e\Gj-J2e=j PeX e (Gj) 

]g e -[E"=j+i a e]Gj-J27=j+i Pe x e(Gj) 


Rj> 0 or Xj cpd Ps 


X 


n— 1 


n 


E e _ [E"=j a e]Gj-J2e=j Pe x e(Gj~) 

E e -[E?=j+i °‘t]Gj-'Le=j+i fa x e(Gj-) 


x Ee 


CX-nGn ftnX r , 


Rj= 0, Xj not cpd Ps 


Proof. First observe that the assumptions imply that the terms Xg(Gj) and Xi(Gj—) in the 
formula are nonnegative for t > j, which legitimates the use of the Laplace transforms. Remark 
also that Rj = 0 for i > j whenever Rj = 0, i.e., for some deterministic io we have Ri > 0 for i < io 
and Ri = 0 for i > i$. 

Let us first suppose that Rj > 0 or that Xj is a compound Poisson process. We prove that for 
j = 1 ,... ,n - 1 , 


E e - E" =; / a iGt-YTl = j Pt x t 


Ee - [££=i E"=i Pi x i{Gj) _ 

-r-=-Ee 

E e _ [E?=j+i a e]Gj-J2e=j+i Pe x e(Gj ) 


aeG e -T,? =j+ 


The key observations are that Xj = Xj (Gj ) and that Gi > Gj almost surely for t = j ,..., n 
by Lemma 1,3.11 The fact that Gj is a splitting time by Lemma, 12.1 f i 1 then yields 

Ee - E?=j aiGi-YTi=j Pt x t 

= Ee~t E ^ h x i{Gj) e - Efcj+i atlGt-Gi]- ELi+1 

= Ee - ^"^' _ E e=j Pe x e(Gj) E e ~ Ei=j+i a e[Ge-Gj]-J2e=j +1 A[W-XRGj)] . ( 2 ) 


The latter factor, which is rather complex to analyze directly, can be computed upon choosing 
°tj = /3j = 0 in the above display. 

Repeating this argument for the case Rj = 0 yields with Lemma 12. If ib provided that Xj is not 
a compound Poisson process, 


Eg - E”=j a eGe-J2e=j Pe x t 


Ee-&i “tps-PiXi-Xt-j+i foXtiPi-) 

--t -Ee Efcj+i 

E e _ [S"=ni a e]Gj-'Ee=j+i Pe x e{Gj -) 


a iGi—YTl = j + \ Pi x i 


It is shown in the proof of Theorem VI.5(i) of [I] that Xj = Xj(Gj—) almost surely, and this proves 
the claim. □ 


In the rest of this section, the following assumption is imposed. 
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G For j = 1,..., n — 1, we have 


X j+1 (t) = K j+1 X j (t) + T j+1 (t), (3) 

where (T 2 ,..., T n ) are mutually independent nonnegative subordinators and If 2 - ■ ■ ■, K n are 
strictly positive. 


Note that Assumption G implies X\ -< X 2 X n . Moreover, it entails that for j = 

1,..., n — 1 and I > j, we have 

i 

Xe(t) = K<jX j (t)+ £ 

i=j +1 


where we have set Ifj = uU j+ i I\i and Kj = 1. In other words, Xg can be written as the sum of 
Xj and I — j independent processes, which are all mutually independent and independent of Xj. 

The following reformulation of © in terms of matrices is useful in Section El Let K be the 
upper triangular matrix with element (*, z + 1) equal to K l+ \ for i = 1,..., n— 1, and zero elsewhere. 
Also write Y(f) := (Yi(t),..., T n (t))', where Yi(f) = X\(t). Equation © is then nothing else than 
the identity X(t) = (I — K')~ 1 T(t). The matrix (/ — K , ) -1 is lower triangular, and element (i,j) 
equals K l - for j > i. 

The cumulant of the subordinator T j ( t ) is defined as 

ej(l3) := -loglEe-^^ 1 ) 


for (3 > 0 and j = 2,..., n. 

The following theorem expresses the joint Laplace transform of (X, G ) in terms of its marginal 
distributions and the cumulants 9 T . However, except for trivial cases, the Laplace transform is not 
the product of marginal Laplace transforms. Still, it can be expressed in terms of these marginal 
transforms in a product-type manner. We call this a quasi-product form. 


Theorem 3.1 Suppose that X is an n-dimensional Levy process satisfying D and G. Then for 
any a,/3 € MT, the transform equals 


Ee 


-[£?=; «*+ESLi+i 9j(T,k=e )]<?,-[£"=,■ K fa] x j 


n— 1 

n Ee -[£,%. +1 a^+EJLi+1 «7(EJ?=£^)]Gi-[ESL i+1 Kfa]Xj 


x Ee 


O^nGn 


Proof. Let j be such that Rj > 0 or Xj is compound Poisson. By Assumption G, we then have 
for a € R+, 


Ee -aG,-E,%- Pe x e( G i) = Ee -aG J -[E2= J - Kfa\X j (G j )~Ei =j+ 1 [J2k=e 

= e f e - aG *-fct=i K il 3 t] Xj ( G j)E ' e -EL J+ i[£^^A] T dGH 
= Ee -[“+£L 3+ i oj (Zk= e ^ fc A)]^-[£?=i 




The claim now follows from Proposition 13.II and the fact that Xj(Gj) = Xj almost surely. 

If Rj = 0 but not a compound Poisson process, the same argument gives the joint transform of 
{Xe(Gj—) : l = j,... ,n} and Gj. In the resulting formula, Xj(Gj-) can be replaced by Xj(Gj) 
as outlined in the proof of Theorem VI.5(i) in Bertoin |I]. □ 
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The following corollary shows that Theorem 13.11 not only completely characterizes the law of 
( X , G) under P, but also its law conditioned on one component to stay nonpositive. Indeed, let Pj, 
be the law of {X{G k + t) — X{G k ) : t > 0} for k = 1,n; it can be checked that this measure 
equals Pj, as defined in Section 12.21 in case R k = 0 and S_ k > 0 P-almost surely. Note that Pj. can 
be regarded as the law of X given that X k stays nonpositive. 


Corollary 3.1 For a,/3€ Mil, we have 


, ^ 1 Ee -[EL, + i^+Efc J+2 0j(Er = ^^)]G,' + i-[Efc j+1 ^ +1 A]A' J+1 


j=k 


Ee 


-[ELi+i a <+E"=j+i *7(EIU E? =J -+i Kp e \x j 


Proof. Directly from Theorem 13.11 and ©. 


□ 


In particular, this corollary characterizes the law of the maximum of a Levy process given that 
it stays below a subordinator. It provides further motivation for studying the law of the vector H 
under Pj,. 


4 The Pj,-distribution of H. 

The aim of this section is to find the Laplace transform of the distribution of H under Pjj, under 
the assumption that 0 is a holding point for {X k (t) — X fc (t) : t > 0} under P. 

We try to follow the same train of thoughts that led us to the results in Section [3J This 
analogy leads to Proposition l4.il which does not yet give the Laplace transform of the distribution 
of H under PjL Therefore, we need an auxiliary result, formulated as Lemma 14.11 which relies on 
Appendix ^ Finally, Proposition 14.21 enables us to find the Laplace transform of the distribution 
of H under Pj,. 

As in the previous section, additional assumptions are imposed on the Levy process X. Here, 
they are significantly more restrictive. The following Assumption H plays a similar role in the 
present section as Assumption G in Section 01 Note that it implies X\ -< Xi X n . 

H Let n = {n(t) : t > 0} be a compound Poisson process with positive jumps only. For each 
j = 1,..., n, we have 

Xj(t) = n(t) — Cjt , 

where Cj decreases strictly in j. 

In the remainder of this section, we write A G (0, oo) for the intensity of jumps of n. We also 
set p^ := sup{i?fc } < Rn^} and cr[ n) := sup{Sj^ : < R^}- In particular, p^ = RrP 

and On ' 1 = . Also, we write for f3 > 0 and i = 1,... ,n, 

ipi(P) := logEe-^ 1 ) 

for the Laplace exponent of — AQ. Since we assume D, is strictly increasing on M + , see the proof 
of Corollary VII.2 of Bertoin [I], Therefore, we can define as the inverse of The function <!>,: 
plays an important role in this section. 

Recall that we used n splitting times to arrive at Proposition 13.11 Here, we only know that 
H k is a splitting time for X under IP[, (see Lemma E3J. In general, however, Hi (i < k ) is not a 
splitting time under Pjj,, and the similarity with Proposition 13.1 1 is lost. 
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Proposition 4.1 Suppose the Levy process X satisfies D. For 7 e , we have 

Eie - Xj=i ljHj = --Ee~ 7j ( p *° . 

' A + £ = 1 7, 


Proof. Lemma l2.HI yields 


Ei e -£i=i ljHj = Kie~^= 1 ' r ^ Hk Eie~^=^ j< ' Hj ~ Hk \ 


In the discussion following o. we have seen that there is a simple sample-path correspondence 
between the laws Pi and P. This yields immediately that Hf. is exponentially distributed under Pj, 
with parameter A. It also gives that the Pl-distribution of {Hj — : j = 1,..., k — 1} is the same 

as the P-distribution {p^ — p * l 'p : j = 1,... ,k — 1}. □ 

Motivated by the preceding proposition, we now focus on the calculation of the distribution of 
the p^ — pj k ^ (that is, their joint Laplace transform). For this, we apply results from Appendix^! 

The following lemma is of crucial importance, as it provides a recursion for the transform of 
— pf : j = 1, — 1} and {p'- 1 — a-'" : j = 1 ,...,*} in terms of the transform of the same 

family with superscript (i — 1). The transforms of the marginals p and p\ l * * _^ also 

appear in the expression, but these transforms are known: for 7 > 0 , i = 1,... ,n (cf. the proof of 
Proposition mu, 

( « (*)\ 

AEe A Pi ^ ) = A+ 7 -^( 7 ). (4) 


Lemma 4.1 Suppose that X is an n-dimensional Levy process satisfying D and H. Then for any 
i = 2,..., n, (3 € M l * + -\7 € M+, we have the following recursion: 


Ee- & 0 -XU v 

( (i) (i)A 
fii-i + AEe -7i v Pi ) 




Pi— i + AEe 
x E e -XUPi{f>W-P'i 






(i) (i) 

Proof. Fix some i = 2,... ,n, and consider the process 7Q_i between a\ and p) . There are 

several excursions (at least one) of the process {Xj_i(t) — fit) : t > 0 } away from 0 between 

a and pf\ and we call these excursions the (i — l)-subexcursions. Each (i — l)-subexcursion 
contains excursions of the processes {Xfit) — Xfit) : t > 0} for i < i — 1; we call these the 

^-subexcursions. To each (i — l)-subexcursion, we assign 2i — 4 marks, namely two for each of 
the i — 2 types of further subexcursions. The first mark corresponds to the length of the last l- 
subexcursion in the (i — l)-subexcursion, and the second to the difference between the end of the last 
Asub excursion and the end of the (£ + l)-subexcursion. Observe that these marks are independent 
for every {i — l)-subexcursion between cr 7 ' 1 and pf\ and that their distributions are equal to those 
of {pf ^ — of ^ : £ = 1 ,..., i — 2 } (the first marks) and {pf + i ' 1 — p^ ^ : £ = 1,... ,i — 2} (the 
second marks). 

The idea is to apply Proposition mu to the process 

Z(x) := inf |t > 0 : X*_i ~ x i-l (vf 1 +tj = xj - ^ X _ ^ , 
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Figure 1: Excursions of {JQ_i(i) — : t > 0} correspond to jumps of Z. 


see Figure El In this diagram, excursions of {Xj_i(f) — X,_, (t) : t > 0} correspond to jumps of Z. 
The relevant information on the subexcursions is incorporated into Z as jump marks. 

Observe that Z is a compound Poisson process with negative drift l/cj_i — l/(cj_i — cf) and 
intensity A/cj_i, starting with a (marked) jump at zero. The jumps of Z correspond to ( i — 1)- 
excursions, and the above marks are assigned to the each of the jumps. In terms of Proposition I A. 11 
it remains to observe that — pf\ and pf 1 — af- 1 correspond to (r_ — Tjv_ )/Cj_ i and r_ / (cj_i — Cj) 
respectively. □ 

With the recursion of Lemma f4.ll we can find the joint transform of p^f 1 — p^ for j = 1,, k— 1, 
which is required to work out Proposition 14.11 This is done in m below. It is equivalent to find 
the transform of P^f+\ — pf 1 for j = 1,..., k — 1, which is the content of the next proposition. We 
have also added p^ — for convenience. The resulting formula has some remarkable features 
similar to the formula in Theorem rm Most interestingly, a quasi-product form appears here as 
well. 

For /3 € R^” 1 > 0, and j = 1,..., k — 1, we define 




= Ci 




(A + Pi). 


Proposition 4.2 Suppose that X is an n-dimensional Levy process satisfying D and H. Then for 
any k = 2 ,..., n, /3 G R^T 1 , 7 > 0, we have 


Ee - ft (^1 -r?) 


t— r Pj + AEe [ C ^+ l(/3)+ 7 1 

3 =l Pj + AEe _ l/o (/3)+ bH 
xEe-[ C " (/ 3 )+ ^ 7 l( p ' 1) -^ ) ). 


r O+D) 
T j +1 J 


O') 


Proof. Since for £ = 2,..., i, by definition of C|(/3), 


ci 


- 1 ha + Pi- 1 ) + —chp)= ctm, 


ci 


it follows from Lemma rm that 

Ee- ft (Pi+i-pf) - [ft fc (ft+^r] 


i + AEe C ^ (/3 ’ 7) (^° 


Ee 




7-i)_ o 7- 1 ) 
j 


) [ C t- 1 (ft+-lr 7 ] W-i 1 ’ cr «-i 1) ) /?£_! + AEe 
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The claim follows from this recursion (start with t = k and note that C(?(/3) = 0). □ 

5 Multidimensional Skorokhod problems. 

In the next sections, we apply results of the previous sections to the analysis of fluid networks. Such 
networks are closely related to (multidimensional) Skorokhod reflection problems , which we describe 
first. Subject to certain assumptions, we explicitly solve such a reflection problem in Section l~> . 1 1 
Section 15.21 describes the fluid networks associated to these special Skorokhod problems. 

Let P be a nonnegative matrix with spectral radius strictly smaller than one. To a given cadlag 
function Y with values in M™ such that T(0) = 0, one can associate a cadlag pair (W,L) with the 
following properties (w E K" ): 

51 W(t) = w + Y(t ) + (/ - P')L(t),t > 0, 

52 W(t) > 0,t > 0 and W( 0) = w, 

53 L(0) = 0 and L is nondecreasing, and 

54 £"=i f 0 °° Wj (t) dLj [t) = 0. 

It is known that such a pair exists and that it is unique; see Harrison and Reiman m for the 
continuous case, Robert M. or Whitt 13?., Thm. 14.2.3] for the cadlag case, and Kella |21] for a 
more general result. 

It is said that (IT, L) is the solution to the Skorokhod problem of Y in M" with reflection matrix 
/ — P' and initial condition w. 

In general, the pair (IT, L) cannot be expressed explicitly in terms of the driving process Y. 
with the notable exception of the one-dimensional case. However, if the Skorokhod problem has a 
special structure, this property carries over to a multidimensional setting. 

5.1 A special Skorokhod problem. 

It is the aim of this subsection to solve the Skorokhod problem for the pair (IT, L ) under the 
following assumptions: 

N1 P is strictly upper triangular, 

N2 the j-th column of P contains exactly one strictly positive element for j = 2,..., n, and 
N3 Yj is nondecreasing for j = 2,..., n. 

In Section 15.21 we show that these assumptions impose a ‘tree’ structure on fluid networks. 
Theorem 5.1 Under Nl—N3, the solution to the Skorokhod problem ofY in M" is given by 
L(t ) = 0 V sup [-(I - P'y'Yts) - (I - P')- 1 ™] , 

0 <s<t 

W(t) = w + Y(t) + (/ — P')L(t), 
where the supremum should be interpreted componentwise. 
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Proof. As W is determined by L and SI, we only have to prove the expression for L. By 
Theorem D.3 of Robert EE, we know that Lj satisfies the fixed-point equation 


Li(t) = 0 V sup [(P'L)i(s) - Wi - Tj(s)] 

o <s<t 


(5) 


for i = 1,..., n and t > 0. 

As a consequence of Nl, we have (I — P')^ 1 = I+P'+.. .+P /n_1 , and the j-tli row of (I — P / )~ 1 
is the j-th row of I + P’ + P /2 + ... + P /J . Therefore, the theorem asserts that 


Li(t ) = 0 V sup 

o <s<t 


i —1 

2^ [p' k Y(s) + P’ k w) 

k =0 


( 6 ) 


The proof goes by induction. For i = 1, © is the same equation as ©■ Let us now suppose that 
we know that © holds for i = 1, ... ,j — 1, where j = 2,..., n. Furthermore, let j* < j be such 
that pj*j > 0; it is unique by N2. Equation © shows that 


Lj(t) = 0 V sup [pj*jLj*(s) — Wj — Yj(s)\ 

o <s<t 



( 



0 V sup 

0 V sup - 2^ Pj*j 

P' k Y(u) + P' k w 

i 

3 

1 

0 <s<t 

_ \ °^ s k =0 


3 1 


= 0 V sup 

0<s<t 


i *- 1 

sup — pj*j P ,k Y(u) + P ,k w —Wj—Yj(s ) 

0<M<S “ L ij* 


= 0 V sup sup 

0 <u<tu<s<t 


k =0 


= 0 V sup 

0 <u<t 


\p' k + x Y(u) + P 

k =0 

[p' k Y(u) + P' k w 


w 


J 3 


Wj - Yj{a) 


k =0 


(7) 


( 8 ) 


where we have used N3 for the equalities © and ©. 

The proof is completed after noting that the j'-th row of P' k only contains zeroes for k = 
j* + l,...,j-l. □ 

Instead of working directly with W, it is often convenient to work with a transformed version, 
W := (/ — P')^ 1 W. The process W lies in a cone C, which is a polyhedron and a proper subset of 
the orthant M!L. Under the present assumptions, at least one edge of C is in the interior of M” and 
at least one is an axis. Below we give an interpretation of W. 

We next establish a correspondence between the event that Wj(t) = 0 and Wj (t) = 0 under an 
additional condition. 


Proposition 5.1 Suppose that Nl—N3 hold, but with ‘nondecreasing’ replaced by ‘strictly increas¬ 
ing’ in N3. Then we have Wj(t) = 0 if and only if Wj(t) = 0, for any j = 1,... , n and t > 0. 

Proof. For j = 1 we have Wj(t ) = Wj(t), so the stated is satisfied; suppose therefore that j > 1. 
Since the matrix (/ — P')^ 1 is lower triangular and nonnegative, we straightforwardly get that 
Wj(t) = 0 implies Wj(t) = 0. 
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For the converse, observe that under Nl—N2 (see the proof of Theorem 15.II we use the same 
notation) 

j-i _ „ j* 


W i (t) = J2[P' k w}.(t) = J2[ 


k=0 


k=0 


I P' k W (t). 
-If 


An induction argument shows that it suffices to prove that Wj(t) = 0 implies Wj*(t) = 0. To see 
that this holds, we observe that by SI and Wj(t) = 0 is equivalent to 


Pj*jLj*(t) — Wj — Yj(t) = 0 V sup [pj*jLj*(s) — Wj — Y)(s)]. 

o <s<t 

The right-hand side of this equality is clearly nondecreasing. Therefore, since Yj is strictly increasing 
by assumption, we conclude that dLj*(t ) > 0, which immediately yields Wj*(t ) = 0 by S4. This 
completes the proof. □ 


5.2 Levy-driven tree fluid networks. 

In this subsection, we define a class of Levy-driven fluid networks, which we call tree fluid networks. 
We are interested in the steady-state behavior of such networks. 

Consider n (infinite-buffer) fluid queues, with external input to queue j in the time interval 
[0, t] given by Jj{t). We assume that J = {J(t) : t > 0} = {(Ji(t),..., J n (t )) : t > 0} is a cadlag 
Levy process starting in J(0) = 0 € M”. The buffers are continuously drained at a constant rate 
as long as there is content in the buffer. These drain rates are given by a vector r; for buffer j, the 
rate is rj > 0 . 

The interaction between the queues is modeled as follows. A fraction pij of the output of station 
i is immediately transferred to station j, while a fraction 1 — Ylj^iPij leaves the system. We set 
Pa = 0 for all i, and suppose that YhjPij — T The matrix P = {pij : i,j = 1,... ,n} is called 
the routing matrix. We assume that for any station i, there is at most one station feeding buffer i, 
and that p^ = 0 for j < i. The resulting network can be represented by a (directed) tree. Indeed, 
the stations then correspond to nodes, and there is a vertex from station i and j if p t j > 0. This 
motivates the name ‘tree fluid networks’. We represent such a fluid network by the triplet (J, r, P ). 
Note that P satisfies Nl—N2 by definition of a tree fluid network. 

The buffer content process W and regulator L associated to the fluid network (J, r, P ) are defined 
as the solution of the Skorokhod problem of 

Y(t) := J{t) - (I - P')rt 

with reflection matrix I — P'. The buffer content is sometimes called the workload, explaining the 
notation W. Importantly, the dynamics of the network are given by SI—S4, as the reader may 
verify. The process Lj can be interpreted as the cumulative unused capacity in station j. 

Associated to the processes W and L, one can also define the process of the age of the busy 
period : for j = 1,... ,n, we set 


Bj(t) := t — sup{s < t : Wj(s) = 0}, (9) 

and let B(t) = ( B\(t ),..., B n {t))'. Hence, if there is work in queue j at time t (that is, Wj(t ) > 0), 
Bj(t) is the time that elapsed after the last time that the j-th queue was empty. If there is no work 
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in queue i at time t, then B; t (t) = 0. Similarly, one can also define the age of the idle period for 
j = l,...,n: 

Ij(t ) := t — sup{s < t : Wj{s) f 0}, 

and the corresponding vector I{t). As a result of these definitions, Ij{t) > 0 implies Bj(t) = 0 and 
Bj(t) > 0 implies Ij(t) = 0 for j = 1,..., n. The quantities Bj(t) and Ij(t) are defined similarly, 
but with Wj replaced by the j-th element of W = (I — P')~ l W. 

The random variables Wj, Bj, and Ij have a natural interpretation. Indeed, let us consider all 
stations on a path from the root of the tree to station j. The total content of the buffers along this 
path is then given by Wj. Consequently, Bj and Ij correspond to the ages of the busy and idle 
periods of this aggregate buffer. 

In the rest of the paper, we assume that the tree fluid network has the following additional 
properties: 

T1 If pij > 0, then p tj > rj/r i: 

T2 Jj(t ) are nondecreasing for j = 2,..., n, 

T3 J is an n-dimensional Levy process, and 
T4 J is integrable and (I — P / ) _1 EJ(1) < r. 

An important consequence of T1 and T2 is that Y is componentwise nondecreasing, except for Y\. 
Consequently, if T1 and T2 hold for a tree fluid network, then Nl—N3 are automatically satisfied 
for the associated Skorokhod problem. Hence, Theorem 15.1 1 gives an explicit description of the 
buffer contents in the network. Note that T4 ensures stability of the network. 

Let us now define the process 

X(t) := (/ - P , ) _1 y(t) = (I - P’)~ l J(t) - rt. 

In view of assumption Tl, the down-stream buffer contents always grow when one of the up¬ 
stream buffers is nonempty. Moreover, under Tl, W is itself a reflected process, that is (W,L) 
is the solution to the Skorokhod problem for X with reflection matrix I and initial condition 
(I — P')~ 1 w. Therefore, each coordinate of IT is a one-dimensional reflected process. A similar 
assumption facilitates the analysis in ffllh Tl . Thm. 4.1 and Lem. 4.2], and [25] . 

In the next proposition, we find the steady-state behavior of the buffer content and the age of 
the busy (and idle) period for the Levy-driven tree fluid network ( J,r,P ). We also consider the 
case where the inequality p^ > Vj/ri in Tl holds only weakly (i.e. pij > rj/ri), as this plays a role 
in priority fluid systems (see Section ROl belowl. 

Recall the definitions of G = G x and H = H x in Sections IP and IP respectively. 

Proposition 5.2 Suppose that T1—T4 hold for the tree fluid network ( J,r,P ). 

(i) For any initial condition IT(0) = w, the triplet of vectors (W(t), B(t), I(t)) converges in 
distribution to ((/ — P')X,G x , H x ) as t —» oo. 

(ii) If the second inequality in Tl holds only weakly, then for any initial condition W(0) = w, 
the triplet of vectors (W(t ), B(t ), I (f)) converges in distribution to ((/ — P')X,G X , H x ) as 
t —> oo. 
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Proof. Throughout this proof, a system of equations like m is abbreviated by B(t) = t — sup{s < 
t : W(s) = 0}. 

We start with the proof of (ii). By Theorem 15.11 we have for any t > 0 

W(t) = [x + A"(f)] V sup [X(t) — A(s)], 

o <s<t 

where x = (I — P')~ 1 w. Moreover, as a consequence of Proposition l5.il we have 

B{t) = t — sup{s < t : W(s) = 0} 

= t — sup |s<t:x + X(s)=0A inf [x + X(u)] 

[ 0 <u<s 

= t — sup <s<i:x + A(s)=0A inf [x + X (u)] ^ , 

^ 0 <u<t J 

where the last equality is best understood by sketching a sample path of X. The supremum over 
an empty set should be interpreted as zero. 

This reasoning carries over to idle periods: 


I(t) = t — sup <s<t:x + X(s)^0A inf [x + A(«)] 

^ 0<u<s 


Due to the stationarity of the increments of {X(t), t > 0} (T3), we may extend X to the 
two-sided process {X(t), t G R}. This leads to 


/ W(t) \ 

( 

B(t) 

=d 

V m y 

V 


[x - X(-t)] V sup - t < s < 0 [~X(s)] 

-sup {s:-t<s< 0 ,-X(s) = [x- X(-t)} V sup_ t < u < 0 [-X(u)]} 


Since x — X(—t) —> —oo almost surely by T4, this tends to 


( sup s < 0 [-X(s)] \ 

- sup { s < 0 : - X (s) = sup u < 0 [-A(u)]} , 

\ -sup{s < 0 : -x(s) + sup„< s [-A(u)]} J 


a vector that is almost surely finite, again by T4. By time-reversibility (see Lemma II.2 of 
Bertoin |3]), the latter vector is equal in distribution to (X,G X , H x ). 

The first claim follows from (ii) after noting that B{t) = B[t) and I(t ) = I(t) by Proposition 

rm □ 

We remark that the above proof does not use T3 to the fullest. Indeed, for the proposition to 
hold, it suffices that J has stationary increments and that it is time-reversible. 

Let us now suppose that the initial buffer content w is random. Proposition 15.21 shows, after a 
standard argument, that {W(t)} is a stationary process if W(0) = w is distributed as //*, where n* 
is the distribution of (/ — P')X. We now show that this stationary distribution is unique. 


Corollary 5.1 Suppose that Tl—T4 hold for the tree fluid network (J,r,P). Then fi* is the only 
stationary distribution. 
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Proof. Suppose there exists another stationary distribution p,Q ^ p*. Let VL 0 * be the corresponding 
stationary process. For any Borel set B in M” and any t > 0, we then have P(VFq(0) £ B) = 
P(W 0 *(i) £ B). Therefore, 

P(W 0 * (0) £ B) = lirn P(W 0 *(t) £ B) 

t—> OO 

/*oo 

= lirn / P(W 0 *(t) £ B\W£(0) = tc)P(IF 0 *(0) £ dw) 

00 Jo 

lim P(W 0 *(t) £ B|W„*(0) = w)P(W o *(0) £ dw) 

t —»oo 

P((J - P')X £ B)P(VF 0 *(0) £ dw) = P((/ - P')X e s), 

where the second last equation is due to Proposition 15.21 This is clearly a contradiction. □ 

Corollary 15.11 answers, for the special case of tree fluid networks, a question from the paper of 
Konstantopolous, Last and Lin EH! on the uniqueness of the stationary distribution. Note that 
for the queueing problem related to ( J,r,P ), the uniqueness of the stationary distribution was 
discussed in Kella E2- In contrast to the setting in ED, we allow for the first component of J{t ) 
to be a general Levy process. 

In the next section, we combine Proposition 15.21 with the results given in Sections 01 and 01 to 
study particular networks. 

6 Tandem networks and priority systems. 

In this section, we analyze n fluid queues in tandem, which is a tree fluid network with a special 
structure. We also analyze a closely related priority system. 

The tandem structure is specified by the form of the routing matrix: we suppose that P is such 
that Pi t i +1 > 0 for i = 1,..., n — 1, and p VJ = 0 otherwise. Observe that we allow Pi ,%+1 > 1, and 
that it is not really a restriction to exclude Pi.i+i = 0; otherwise the queueing system splits into 
independent tandem networks. 

In all of our results, we suppose that the tandem system ( J,r,P ) satisfies Tl—T4. We rule 
out the degenerate case where the first j > 1 components of J are deterministic drifts, since an 
equivalent problem can then be studied with the first j stations removed. We also impose the 
following assumptions on the input Levy process J: 

T5 J has mutually independent components, and 

T6 The Levy measure of J\ is supported on M + . 

Observe that under T2—T3, T5 implies that J2 ,..., J n are independent nonnegative subordi- 
nators. 

This section consists of three parts. In Section RT~T1 we are interested in the joint (steady-state) 
distribution of the buffer contents and the ages of the busy periods for fluid tandem networks, i.e., 
in the distribution of (PF(oc), B(oo)). Section RT~?1 considers the situation of a single compound 
Poisson input to the system. For that system, we are also interested in the ages of the idle periods, 
i.e., in the vector /(00). In Section lb.,‘II we analyze buffer contents and busy periods in a priority 
system. 
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6.1 Generalities. 

To find the joint distribution of VF(oo) and B{ oo), throughout this section denoted by W and B 
respectively, we rely on Proposition QT21 This motivates the analysis of X(t) = (/ — P')~ 1 J(t) —rt. 
For i = 2,... ,ra, we dehne the cumulant of Jflt ) by 9](ft) := — logEe - ^ 3,7 ^ 1 ), (3 > 0. As in Section^ 
we write ifi (defined by = logEe - ^ 5 ^ 1 )) for the Laplace exponent of —Xj. Its inverse is again 

denoted by <!>,. 

Under T2 and T6, the Levy measure of X is supported on Rif. Moreover, as we ruled out trivial 
queues in the network, each of the components of X has a nondegenerate distribution. Therefore, 
let us recall that the following holds (see, e.g., Theorem VII.4 in 0): for at,/3> 0, (a,/?) ± (0,0), 
f3 <3?j(a), i = 1,..., n, we have 

E -aGi-m = _E.Y f (l) $ * ^ . ( 10 ) 

This identity plays a crucial role in the results of this section. For notational convenience, we shall 
write that cud holds for any a, (3 > 0 , without the requirements (a, (3) (0,0) and f3 4>j(a). 

Now we can formulate the main result of this subsection. We remark that the first formula also 
holds if J\ is not necessarily spectrally positive. For instance, it allows for phase-type downward 
jumps; see d for the joint transform of Xj and Gj in that case. 


Theorem 6.1 Consider a tandem fluid network ( J,r,P ) for which T1—T6 holds. Then for uj, (3 £ 
Rif, the transformequals 


]E e -[E?=j+i ( u e)+J27=j+i(pe-ure-i—re)ui£+Yii = j Pe\Gj—u>jXj 
Ee _ [^<=j'+i e e^ u} B+'Ee=j+i(Pe-ie r e-i~ r e)^e+Yi , e=j+i Pt\Gj~Pj,j+iVj+iXj 


n— 1 

n 


x Ee 


dnGn LdnXrj 


Consequently, we have for w,f3 £ 




E = -EV n (l)^ (/?n) 7 ^ 

Pn ^Pn\^n) 

t} (E”=j+1 ( u i) + E?=j+i(Pe-iJ r e-i ~ r t) UJ i + E™ =J 


>< n 

3 = 1 


U (E"=j+i 6 1 M + E"=j+i(p<-i,< r <-i - n)wi + E"=j+i Pe) - Pjj+iUj+i 

n_1 E"=,-+1 e i (vt) + E"=,-+i(P*-M r *-i - r e )we + E"=i+1 Pt ~ ^(PjJ+iUj+i) 


E”=i+i e i ( w <) + E"=j+i(p*-i,< r <-i - n)^t + ELi Pt - 

Proof. By Proposition I5.2f ib (W,B) =<j ((/ — P')X,G x ). Hence we have 

Ee -(u;,W)-(f3,B) = e e -({I-P)u,i!-P , )- 1 W)-(fi,B) = Ee -(/3,G)-((/-P)o;,X>_ 


( 11 ) 


Now note that the stability condition T4 for (J, r, P) implies D for X by the law of large 
numbers. Thus, in order to apply Theorem 13.11 for CD, it is enough to check that G holds. 
Standard algebraic manipulations give 


Xi(t) = Ji(t) - rfl 


and 

X i+ i(t) = Pi,i+iXi(t) + Ji+i(t) + (Pij+m - r i+ i)t 
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for i = 1, ..., n — 1. Hence, G holds with Ki = and Yj(t) = Jflt) + — rf)t. 

As a result, we know that from Theorem cm 


Ee -</3,G>-((/-P)w,x> = Ee -</3,G)-<S,x> 

^ Ee -[E"=j+i ej (Ek=e «JS3fc)+E” =i A]Gi-(EZ=i 
Ee -[E? =3+ i(£*=< K t^)+E?=,+i A]G J -(EZ= i+1 


x Ee 


/5n Gn LOn -V ,j 


where we have set uj = (I — P)w for notational convenience. 

The reader may check that Ylk=j Kj^k = Vj and Ylk=j+i Kj '= Pj,j+iWj+i, leading to the 
hrst claim. The second assertion is a consequence of the first and (II01) . □ 

Theorem 16.11 extends several results from the literature on the steady-state distribution of the 
buffer content for tandem Levy networks. In particular, if J(t ) = 0)', P = (pij), with pi 2 = 1 

and zeroes elsewhere, if one chooses (3\ = = 0 and lo\ = 0 in Theorem 16.11 then one obtains 

Theorem 3.2 of D§bicki, Mandjes and van Uitert [!J|. Additionally, if one chooses fi\ = /f ?2 = 0 and 
supposes that ,J\ is a subordinator, we recover the results of Kella E03- 

Even if the Laplace transform of ( Gj,Xj ) can be inverted, it is generally not straightforward to 
invert the Laplace transform of (W, B ) given in Theorem 16.1 1 Some progress has been recently made 
in case n = 2; for a Brownian fluid system, Lieshout and Mandjes m calculate the distribution 
of W. Avram, Palmowski, and Pistorius [2] study a compound Poisson setting with exponential 
jumps. A different type of explicit solution can be found in the work of Harrison M; he gives an 
example closely related to the framework of the present paper. 

For use in Section f6.3l we point out that the expression in Theorem 16.1 1 is Ee B ) if the 

second inequality in T1 is weak, cf. Proposition 15. 2f iih 


The lengths of the busy periods. 

Besides the Laplace transforms of the ages B of the busy periods, Theorem 16.1 1 also enables us 
to find the Laplace transforms of the length V of the steady-state running busy periods. Indeed, 
let Di, i = 1 ,... ,n denote the steady-state remaining lengths of the running busy period, so that 
Vi = B, + D, ■ We know that I), and B, are equal in distribution. In fact, following for instance JTJ 
Sec. V.3], we have 

(Bi, Di) = d (UiVi, (1 — Ui)Vi), (12) 

where Ui are i.i.d. and uniform on [0,1]. 

For the Brownian (single-station) fluid queue, the following result is Corollary 3.8 of Salminen 
and Norros |36| . 


Corollary 6.1 

0, a + (3, 


Consider a tandem fluid network (J, r, P ) for which Tl—T6 holds. 


Ee -aBi-PDi 


—EXj(l) 


$j(oQ -$iQ 3 ) 

a — /3 


Then for a, fl > 


Moreover, we have for a >0, 


Ee~ aVi = 


—EXj(l) 


d$i(a) 

da 
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Proof. Since the second claim follows straightforwardly from the first, we only prove the first 
expression. Following m, we have for a f3, 


{a - P)Ee- aBi ~P Di 


(<a - p)Ee^ a -^ UiVi ^ Vi = ( a - @)E [ e-^ a ~^ uVi ^ Vi du 

Jo 

e[ e~ uVi du = e[ e~ uVi du-E [ e~ uVi du. 

J/3 Jo Jo 


The two identities that result upon setting j3 = 0 and a = 0 can be used to express the first and 
second expectation in terms of the Laplace transform of Bi and Di respectively; this yields for 
a ^ /3 


E e ~ aB i-P D i = 


a — (3 


aEe 


—olB' 


1 - f3Ee~ 0Bi 


where we have used the equality in distribution of Bi and Di. Application of m completes the 
proof. □ 


6.2 A single compound Poisson input. 

In this subsection, we examine a tandem fluid network with a single compound Poisson input j25| . 
The following assumption formalizes our framework. 


T7 Pi : i +1 = 1 for i = 1,,n — 1, while Pij = 0 otherwise, and 

T8 Ji is a compound Poisson process with positive drift d and intensity A, and Jj = 0 for j = 
2,... , n. Moreover, rj decreases strictly in j and EJ(1) < r n . 

An important consequence of T7 and T8 is that 

(rj - r k )uj = - ipk(v), (13) 

which simplifies the resulting expressions in view of fact that we often deal with ratios of the 
fluctuation identity m- Interestingly, it is also possible to study (joint distributions of) idle 
periods under these assumptions. 

The following corollary collects some results that follow from T7 and T8 and Theorem 16.11 
Many interesting formulas can be derived, but we have selected two examples for which the formulas 
are especially appealing. 


Corollary 6.2 Consider a tandem fluid network ( J,r,P ) for which T7—T8 holds. 


(i) For i = 1,..., n, and u>, ft > 0, we have 




—EXj(l) 


®i{(3) ~ u 
(3 + (r;_i - n)u 


$i_i((r-j_i - n)u + P) 
$i_i ((ri_i - n)uj + (3) - u ' 


Moreover, P(Wi = 0) = P {Bi = 0) = 
(ii) For i = 2,..., n, u, (3 > 0, we have 


E 


e-uWi-m. Wil 



EXj(l) 

d - n -1 $j_i((rj_i - n)u + (3)-u' 
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Proof. To prove (i), apply Theorem 1(1.II to obtain for i = 1,,n, 


E = Ec _ 8g ,_„ x . 

Ee ~ [(ri_i -ri)u+l3\Gi-i -uXi -1 

With (II01) . this leads immediately to the given formula after invoking (II.'ll) . 

We find P(W* = 0) upon choosing u> = 0 and noting that 

P (Wi = 0) = P (Bi = 0) = lirn Ee _|3Gi = -EXfll) lim = EXi ^\ 

/3—kx> /3—>oo f3 d — Vi 

where the last equality follows from Proposition 1.2 in :3|. 

The second claim uses a similar argument; it follows from Theorem 16.11 that for i = 2,. .. ,n 

Ee -w i W i -/3 i _i.B i _ 1 -/3 i B i = _— E(! 3iGi 

Pg—[(r»-i—ri)a;i+/3j]Gi_i—cjjXj-i 

and the numerator of the fraction tends to P(Wj_i = 0 ) as /%_i —► oo. Now apply (fTTTl) and m- 
□ 

We end this subsection with an application of the theory in Section 0 which enables us to study 


the idle periods in a tandem fluid network satisfying T7—T8. For 7 £ 


■ok — 1 


we set 


^•(7) := Cj 


s Pr-aB”) 


which is similar to the definition of Cj in Section 

Proposition 6.1 Consider a tandem fluid network (J, r, P) for which T7—T8 holds. For 7 £ M™ , 
we have 

n 

Ee -< 7 ,/> = 1-J2 ¥ (. w k = 0)Ej|; e“£*=i 11He (1 - e~ lkHk ) 
k=1 

where P (Wj = 0) is given in Corollarv \6.2tf i). and 


E [ :e -Eti7^£ = 




A + ELi7< 

X g A + Efci 1 It - Efc /+1 > - c fc $ i+1 (PjVi(7)) 


(14) 


A + Eti7^-Et J+ it7t-c^(^(7)) 

Proof. Note that T7 and T8 imply H. The first claim follows from Proposition 15.21 and the facts 
that for k = 2 ,..., n, 

Ee - Ef=i7 eH e = Ee-Efei + E i (i - 7 fe f7 fc )l P(Z fc = 0), 

and Ee 71 ^ 1 = 1 — e[ [l — e -71 ^ 1 ] P(A"i = 0). These identities follow after observing that H^ 
vanishes on the event { X = 0 }, and that { X^ = 0 } is the complement of { X & > 0 }. 

Let us now prove the expression for the Pjpdistribution of (Hi,..., Hf)'. From Proposition l4.il 
and Proposition 14.21 we know that 


E l e -ELi7^ = 


AEe 


^(7)(pi 1} <T ^ 1) ) ^ Y3 e=1 'ye + AIEe ^+ + i 1) ) 


n 


A + E*=i7 1 j=1 J2i = i 7 1 + AEe 
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The proof is finished after invoking 0 and noting that for j = 1,..., k — 1, 




x +^2,it+ vk j( 7 ) 


*=i 


Cfc 

Cj + 1 


A + J^7£ + ^+i(7) 


€=1 


k —1 /c—1 

A + X^ _ X] “7/. 


t =1 £=j+i 


Q 


and 


£fc 
Cl L 


fc —1 / c —1 

A + V\ (7)] = A + ^ 7£ - ^ —■ ~u, 

J Co 

1=1 t=\ *■ 


as the reader readily verifies. 


□ 


6.3 A priority fluid system. 

In this subsection, we analyze a single station which is drained at a constant rate r > 0. It is fed by 
n external inputs (‘traffic classes’) J\(t),... ,J n (t), each equipped with its own (infinite-capacity) 
buffer. The queue discipline is (preemptive resume) priority , meaning that for each i = 1 ,n, 
the i-tli buffer is continuously drained only if first i — 1 buffers do not require the full capacity r. 
We call such a system a priority fluid system. 

The aim of this section is to find the Laplace transform of ( W,E ), where Wj = Wj(oo) is the 
stationary buffer content of class-j input traffic, and Ej = Ej( 00 ) is the stationary age of the busy 
period for class j. We impose the following assumptions. 

PI J is an n-dimensional Levy process with mutually independent components, and its Levy 
measure is supported on M”, J(0) = 0, 

P2 Jj(t) are nondecreasing for j = 2,..., n, and 

P3 J is integrable and )T” =] E,/,;(l) < r. 

The central idea is that IT evolves in the same manner as the solution to the Skorokhod problem 
that corresponds to a tandem fluid network (J, r, P), with r = (r,..., r)' and P = ( pij) such that 
Pi,i+i = 1 f° r * = 1, ■ ■ ■ , n — 1 and p^ = 0 otherwise. This equivalence has been noticed, for instance, 
by Elwalid and Mitra m- It allows us to use the notation of Section 16.11 

It is important to observe that PI—P3 for the priority system implies Tl—T6 for the cor¬ 
responding tandem fluid network, except that the second inequality in Tl only holds as a weak 
inequality. However, as remarked in Section r6.ll the Laplace transform of the distribution of (W, B ) 
is then still given in Theorem 16.II 

The steady-state ages of the busy periods E can also be expressed in terms of the solution (W, L) 
to this Skorokhod problem, but it does not always equal B as in Section 16.11 To see this, notice 
that if class-1 traffic (highest priority) arrives to an empty system at time t, we have W 2 (t) = 0, 
while W r 2 (t) > 0 so that B‘>{t) > 0. However, it must hold that ^(f) = 0. 

Still, the following theorem shows that it is possible to express the distribution of (IT, E) in 
terms of (IT, B). 

Theorem 6.2 Consider a priority fluid network for which PI—P3 holds. Then for uj, /3 € M" , the 
transform equals 

n 

Ee -(UJ,W)-(P,B) + JE [e“ PiB t ^ . Wj==Wn = 0 . 

3 =2 
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Proof. In principle, Ej equals Bj, except when Wj = 0. In fact, it follows from the above reasoning 
that 


Ee -(u,W)-(i3,E) = E W2 = =Wn = 0 

n 

+ E [ e ~ E * =1 UiWi ~ E3t=1 Wj > 0, W j+ 1 = ... = W n = 0 

3= 2 


Now use the fact that {Wj > 0} is the complement of {Wj = 0} and rearrange terms. 


□ 


If the J 2 , • • •, J n are strictly increasing, it can be seen (for instance with Theorem 1(1.1 1) that 

'■J - 1 


E 


( 1 - e~h B i ) ; Wj = ... = W n = 0 


= 0. 


Therefore, in that case, we have the equality in distribution (W,E) =d (W,B). 

Another important special case is when Jj,.... ,J n are compound Poisson processes, say with 
intensities Ai,..., X n respectively. Much is known about the resulting priority system, see for 
instance Jaiswal nsi for this and related models. To our knowledge, the distribution of (W, E) has 
not been investigated. However, it is given by Theorem 1(121 and Theorem irm upon noting that 
0f(uj) —> Xi as u —> 00 . Since it is not so instructive to write out the resulting formulas, we leave 
this to the reader. 
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A Appendix: some calculations for a compound Poisson process 
with negative drift 

In this appendix, we study a compound Poisson process Z with negative drift, and derive some 
results on the excursions of Z — Z_ from 0, just before its entrance to 0. These results are applied 
in Sectional 

Let us first fix the notation. Throughout this appendix, Z is a Levy process on (P,^ 7 , P) with 
Laplace exponent 

V>_z(/3) := logEe - ^ 1 ) = c fj - \ I (\ - e ~P z \ F(dz), 

J R+ ^ ' 

where c > 0, A € (0, 00 ), and F is a probability distribution on (0, 00 ). That is, Z is a compound 
Poisson process under P with rate A and negative drift —c, and its (positive) jumps are governed 
by F. We suppose that EZ(1) < 0, so that Z drifts to — 00 . In analogy to Section 0 the inverse of 
tp-z is denoted by it is uniquely defined since ip-z is increasing. Observe that <h_z(0) = 0. 

Set To = 0, and let T* denote the epoch of the i-th jump of Z. To the i-th jump of Z, we associate 
a vector of marks, denoted by M % E (for some m E Z+). We suppose that Mi is independent 
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of the process T = {T n : n > 1}, and that it is also independent of ( Z(Tj) — Z(Tj—), Mj) for j ^ i. 
However, we allow for a dependency between M % and Z(Ti ) — Z(Ti~). In fact, an interesting choice 
for Mi is Mi = Z(Ti) — Z(Ti~) (so that m = 1). 

Define t_ as the first hitting time of zero, and IV_ as the index of the last jump before r_, i.e., 

t_ := inf{i > 0 : Z(t) = 0}, IV_ = inf{n > 0 : Z(T n+ i —) < 0}. 

Write P^ for the law of Z + £ under P with initial mark Mq = M. We suppose that the initial 
condition (£, M) is independent of Z , and has the same distribution as ( Z(T \) — Z(T\—),Mi). 
Observe that both r_ and IV_ are P^-almost surely hnite, and that (by the Markov property) the 
‘overshoot of the first excursion’ T ^_+1 — t_ has an exponential distribution with parameter A. 

In this appendix, it is our aim to characterize the P^-distribution of r_ (excursion length), 
r_ — Tn_ (excursion ‘undershoot’), and M (mark of the last jump). Overshoots and undershoots 
have been studied extensively in the literature. However, as opposed to what we have here, these 
results are all related to the situation that a Levy process can cross a boundary by jumping over 
it (strictly speaking, this is the only case where the terms ‘overshoot’ and ‘undershoot’ seem to be 
appropriate). See Doney and Kyprianou m for a recent contribution and for references. 

In view of the results of Dufresne and Gerber nu, it is tempting to believe that r_ — Tjyr_ has 
an exponential distribution. However, it turns out that this ‘undershoot’ has a completely different 
distribution. 


Proposition A.l We have for f3, 7 > 0 and k € M™, 

„ TV > \P + 7 - C<h-z(7) + A] 

? (3 + 

\J3 + AE ? e-T' r -] Ee-(^ +A )«/ c -< K ’ M > 

(3 + AEe-(^+ A )£/ c ' 

To prove this proposition, we need an auxiliary result on Poisson processes. Consider a Poisson 
point process N(t) with parameter //, and let ( be a positive random variable, independent of N. 
Let Aft) be the backward recurrence time process defined by N, that is the time from £ to the 
nearest point to the left. The following lemma characterizes the joint distribution of IV(£), A(£), 
and C- 


Lemma A.l We have for (3, 7 > 0 and 0 < s < 1. 

E s Ar (C) e -/3^(C)-7C _ (L £ c -(/3+7+m)C _|_ S V‘ £ c -l7+(l-a)^lC j 

(3 + s/i (3 + s/j, 

Proof. We only prove the claim for 7 = 0; the general case follows by replacing the distribution 
of C by the (defective) distribution of ( given by Ee - ^ = Ee“^ +7 ^. Let Ho = 0 and U\, U 2 ,... be 
the location of consecutive points of N. Observe that 


Ea N(C) e -PA(0 


Y 0 < C ~U n < U n+ 1 - U n 

n =0 

00 roc rt 

Y * n / e-^)^P [/n (dx)P c (dt) = Y +P), (15) 

—" J ° J ° n=0 


n =0 


24 








where 


MP) :=E[e-frH(>Un 

Clearly, po(P) = Ee - ' 5 ’’. If we let B be the forward recurrence time process, we have for n > 1, 




-E 


-P(C~u n ). [/ n _i < ( < U n 


= E 
= E 


e -P(C-U- n _i)+/3(l7 n -l7 n _i) ; ^ > [/n _ i 


-E 


-P(G-u„). u n l < ( < U n 


J3(U n —U„-i) 


E 


3 -/3(C-t/ n -i) ;C > Un _ x 


-E 
li- (5 


J3B( C) 


IV(C) = n- 1 P (IV(C) = n-l) 


,_i(/3)-P(IV(C) = n-l)], 


where we used the lack-of-memory property of the exponential distribution for the last equality. 
After iteration, we obtain 


MP) = 




v-P 


n— 1 


Ee^ c - Y 


i =0 




M-/5 


p(iv(c) =»). 


Therefore, taking 0 < s < (3/^ (later we may use an analytic-continuation argument), we deduce 
from m that 


E 


S N( 0 e -/3A(0 


) n oo n—1 


n=0 


P(iV(C) = i). 


The double sum in this expression can be rewritten as 


Sfl 


P + Sfl 


LXJ 

Y « ip (^(C) = *) = -t^-Ee-^ 1 -^, 

jj Sfl 


i =0 


and the claim follows. 

Lemma mu is the main ingredient to prove Proposition mu 

Proof of Proposition IA.1L The crucial yet simple observation is that 

-t n _ )-7 t--{k,m n _ ) 


□ 


E^e 

= Et 


3 -/3(r_-T JV _)-7r_-(K,Mjv_>. jy_ _ q 


+ Et 




= E e _ ^ A+/3+7 ^/ c_ ^ K,M ^ + E^ 


/ 9 ( t _ -T n _)-"/t- - ( k,Mn_ ) 


; AT > l 


(16) 


To analyze the second term, we exploit the fact that there are several excursions of Z — Z from 0. 
Therefore, we set 

C(f) := inf{s > 0 : Z(s) — Z{ 0) = —t}, 

where an infimum over an empty set should be interpreted as infinity. 

It is obvious that C is a subordinator with drift 1/c, and that it jumps at rate A/c with jumps 
distributed as r_ under P^. This observation implies with Theorem VII. 1 of Bertoin |5] that 


$-z( 7 ) = - + -(l-E 5 e- 7T -). 

cc 


( 17 ) 
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Lemma mu can be applied to the Poisson process N constituted by the jump epochs of C, 
fi = X/c, and £ = £. Each jump of C corresponds to an excursion of Z — Z_ from 0, for which the 
‘excursion overshoot’, the excursion length, and the marks of the last jump are of interest. Observe 
that these quantities have the same distribution as r_ — T/v_, r_, and Mjy_ respectively. Using the 
notation of Lemma mu this yields 


Ef 




= E 


(E ? e“ 7T -) N(e) 1 e-PMO/c-rt/c. jv(f) > lj ^ e -P^_^T N _)- 1 r--(K,M N _)_ ( 18 ) 


Therefore, Lemma mu yields 


E 


s A r (U-l e -/3A(0/c-7€/c. jv(£) > i 


E [ s JV(S) e -/3A($)/c-7f/cj _ JE e -0+/3+7)f/c 


E e _ ^ 1_s ) A+7 ^/ c - Ee - ^" 1 "^" 1 " 7 ^/ 0 


As + (3 - 

Upon combining this with csD and CHI , we arrive at 

\/3 + AE c e“ 7T -] E e -(^+/3+ 7 )€/c-<«,M> 


E^e 


—0{t~ —Tff_ )—7T_ — ) _ 


AE ? e-^ T - +p- AEe - W 1_B; e e 7T ")+i)£/c + AEe-( A +^+ 7 )£/ c ’ 

which, with the help of m, reduces to 

\P + 7 - + A] Ee-( /3+7+A ^/ c -< K ’ M > 

(3 + 7 — c^-zin) ~ A (Ee- $ -z( 7 )« - 1) + AEe-^+ 7+A )^/ c ' 

By definition of $>-z, we have 


7 = i/j-z(<&-z('y)) = c$-z( 7) + A ^Ee $ - z(7 ^ _ ) 


and the claim follows. 


□ 
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